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We present a new and exploratory approach to determine the A/9(/9)-shift in the 0(3) nonlinear fr-niodel. The 
method is based on a scaling hypothesis for a free energy difference, which is assumed to be valid in a situation 
where the mass gap correlation length { is of the order or larger than the linear extent L of the considered 
square lattice sizes. The free energy difference arises from the finite volume constraint effective potential of the 
theory. While the constraint effective potential is calculated in numerical simulations employing a variant of the 
multicanonical ensemble on medium sized lattices, it is possible to estimate A/9(/9) up to a value of /? = 2.8 in the 
standard parameterization of the model. 



1. INTRODUCTION 

The D = 2 nonlinear 0(3) sigma model with 
the lattice action 
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is believed to be a prototype model of asympto- 
tically free theories and therefore deserves spe- 
cial interest in view of asymptotically free gauge 
theories. It possesses a finite mass gap and in the 
perturbative regime of couplings /3 its finite mass- 
gap correlation length £,~^ = mcjapa is supposed 
to scale like 

e(/3) = Q^(l + ...). (2) 

Possible corrections to this scaling law (...) have 
been considered in perturbation theory and are 
small. In this paper we will be concerned with 
numerical studies of the cr-model. Over the years 
there have been several pioneering numerical stu- 
dies [1] trying to infer the A/3(/3)-shift of the 
theory. These calculations located the /3-region, 
where perturbative scaling is valid, and attemp- 
ted to determine the mass gap of the theory in 
physical units. Hereby corresponds A/3(/3) to the 
positive shift in (3 implied by a doubling of the 
correlation length and its perturbative prediction 



^«'''> =!!<'+ 27?+ °"/"'»- 



(3) 



Most of the numerical simulations have used mea- 
surements of two-point correlation functions in 
order to determine the mass gap [1] and the qua- 
lity of these simulations was boosted with the in- 
vention of the refiection cluster algorithm. Other 
numerical calculations have resided on the renor- 
malization group approach [2]. Nevertheless the 
exponential increase of ^ with (3 makes direct and 
sensible two point function measurements impos- 
sible for /3-values above /3 Ri 2 on today's compu- 
ters. With the help of the renormalization group 
it was possible to extend the accessible /3-region 
perhaps up to /3 Ri 2.4, but still, very large lattices 
are needed in these simulations. 

The constraint effective potential Uej j {M) 
(CEP), of spin models with continuous global 
0(N) symmetry is a function of the length M 
of the mean-field 
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and is defined by rewriting the finite lattice path 
integral of the theory in the form 



Z = c dMM 
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which can be achieved in principal by introducing 
suitable (5-functions. Note that in the integral a 
term M^~'^ appears. It counts the degeneracy of 
states with given M and is proportional to the 
surface of the #-dimensional sphere. The CEP 
has been extensively studied in another context 



with analytic methods [3] as well with numeri- 
cal simulations [4], namely in _D = 3,4 0(N)- 
symmetric ^'^-theories at a value of the bare quar- 
tic coupling A = cxd. These models are ferroma- 
gnets and possess, unlike the D = 2 0(3) non- 
linear (T-model, a symmetry broken phase with 
a nonzero field expectation value S. In the bro- 
ken phase the CEP has on finite lattices at Mmin 
a minimum, Mmin is finite in the thermodyna- 
mic limit, corresponding to S. Here we may note 
that both analytic as well as the numerical consi- 
derations in these model were concerned with the 
shape of the CEP in the vicinity of its minimum 
and that not much is known about it in regions 
of its argument far away from the minimum, e.g. 
at M = 0, see the discussion below. 




Fig. 1 a plot of the CEP at a value of /3 = 2.4 
on 24^, 36^ and 70^ lattices, as it was determined 
in numerical simulations, which we describe la- 
ter. We have normalized the value of Ueff(M) at 
its clearly visible minimum to 0. But, how then 
knows the CEP of the D = 2 cr-model about the 
symmetric phase of the theory ? Qualitativly we 
make the following observations: Firstly the lo- 
cation of the minimum Mmin moves to smaller 
values if L is increased. Secondly the difference 
^Ueff = Ueff(M = 0) - Ueff(Mmin) decreases 
when L is increased. This second observation is 
essential here, as the difference AUeff plays the 
role of a potential barrier in the CEP for states 
which, at Mmin on the finite lattice are probable 
states, and states, which at M = are suppressed 
in the path integral on finite lattices. These states 
are however characteristic for the theory in the 
thermodynamic limit. After all, there we expect 
a vanishing field expectation value < M >= 0. 
We therefore expect that in this limit 



AUeff = 0, Mmin = 
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for all values of /3 in the D = 2 sigma-model. At 
this point we may come back to above mentioned 
ferromagnets. There it is expected that 



AUeff oc L° 



(7) 



with a = D — 2. Thus the volume normalized 
CEP Ueff/L^ of such theories will in the thermo- 
dynamic limit have a convex shape with constant 
values between and some My,. 



Figure 1. Plot of the constraint effective potential 
at /3 = 2.4 on 24^, 36^ and 70^ lattices. 

The D = 2 0(3) nonlinear cr-model however ex- 
hibits a completely different physical situation. It 
possesses just a single symmetric phase and there- 
fore an appropiately defined field expectation va- 
lue S is identical to zero. Nevertheless, on finite 
lattices, e.g. square lattices with periodic boun- 
dary conditions, one can observe a nonvanishing 
expectation value of the above operator < M > 
even for the sigma-model, when the mass gap cor- 
relation length ^ is of the order or much larger 
than the linear extent L of the considered latti- 
ces. For purpose of illustration we display here in 



2. SCALING HYPOTHESIS 

Considering the path integral of the theory for 
states with M« = (M^, ...,M^) fixed «,e,, intro- 
ducing the constraint partition functions 
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we note that due to 0(#)-symmetry these func- 
tions are functions of M alone. Furthermore a 
free energy difference of two such states can be 
expressed in terms of the CEP in a very simple 
way 
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Ueff{M,)-Ueff{M,). (9) 



Thus AUeff = UeffiM = 0) - UeffiMmin) Cor- 
responds to a free energy difference. In statistical 
mechanics there exists a powerful scaling hypo- 
thesis, Fisher scaling, which is applicable in the 
framework of finite size scaling theory. At the cri- 
tical point I.e., in situations where the correlation 
length ^ is of the order or much larger than any 
linear extent L of the finite volume system, a free 
energy or their difference is a function of the ra- 
tio ^/ L alone. Similar idears have been presented 
in [5] for the case of the Ising model and for the 
helicity modulus of the the D = 2 0(3) nonlinear 
a model [6]. In this paper we will test a corre- 
sponding scaling hypothesis for the D = 2 0(3) 
nonlinear cr-model i.e., we will assume that 



AUeff = AUeff (UL) 
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for values of /3 and lattice sizes L where the con- 
trol parameter ^/L is large. Specifically we will 
assume an analytic form for the /3-dependence of 
^ in the vicinity of a given /3o 



e(/3)=e(/3o)|e^''/^°; 



A/3 = /3 - /3o 



(11) 



and we will consider the coefficient Bq as free pa- 
rameter. Bo = l/(27r) in perturbation theory. Bq 
and also the A/3(/3)-shift are then determined by 
the behavior of AUeff in the vicinity of /3o as 
determined by the Monte Carlo simulation. 

3. NUMERICAL SIMULATIONS 

The CEP can be obtained in numerical simula- 
tions from the probability distribution function 
P(M), which denotes the probabilty of states 
with mean field length M in the canonical en- 
semble 
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However, simulations of the path integeral on lat- 
tice sizes and /3-values considered here strongly 
suppress states with small values of M . We the- 
refore have used a variant of the multicanonical 
ensemble, in which we have added a weight func- 
tion W(M), enhancing the appearance of states 
with small values of M in the simulation. This 
way it was possible to obtain distribution function 



estimates for the function Ueff (M) on the whole 
relevant M-interval. Typically 400 bins were used 
for these distribution functions and Fig. 1 shows 
an example. We have considered four /3- values, 
namely /3 = 1.6, 2.0, 2.4 and /3 = 2.8. The consi- 
dered range of square lattices was between _L = 24 
and _L = 80. We use periodic boundary conditi- 
ons. Each simulation was performed with a stati- 
stics of about 10® sweeps. Error calculation was 
performed with the jackknife method. For each 
pair of values L and /3 we have determined AUeff 
by fits to the distribution functions. The value of 
Ueff at Mmin was determined by fitting a para- 
bola to Ueff in the vicinity of the minimum, in 
the same way a polynomial fit was used to deter- 
mine Ueff(M = 0). 




Figure 2. Plot of AUeff for a 36^ lattice close to 
/3 = 2.4 as a function of /3. The displayed curve 
is a polynomial fit to the data. 

4. RESULTS 

Here we refrain from a detailed discussion of 
our results from simulations at /3 = 1.6,2.0 and 
2.8, the final A/3(/3)-shift for these /3- values can be 
found in Fig. 4. Results at these /3-values appear 
similar to results for /3 = 2.4, which we present 
now. In Fig. 2 we display results for the quantity 
AUeff as a function of /3 close to /3 = 2.4. One 
observes a monotic increase of AUeff with /3. As 
can be seen, the functional dependence on /3 is 



very smooth and therefore a polynomial approxi- 
mation of the form 

AUeff = a + &(/3-/3o) + c(/3-/3o)2; /3o = 2.4(13) 

gives a very good description of the analytic be- 
havior with l3. This analytic form can then be 
employed in order to test the scaling hypothesis 
Eqs. 10, 11. Figure 3 displays the _L-dependence 
of AUeff at /3 = 2.4. Lattice sizes range inbet- 
ween _L = 24 and L = 70. The decrease of AUeff 
with increasing L is observed. The curve in Fig. 
3 corresponds to the rescaled data on 36^ latti- 
ces, Eq. 13 and Fig. 2. We have used the scaling 
hypothesis AUeff = AUeff (^/L) and the ^(/3) 
dependence Eq. 11, with a coefficient Bq as de- 
termined by a x^"fit to the data. We note that 
the scaling hypothesis is well satisfied. Figure 4 
then contains the inferred A/3(/3)-shifts for all our 
/3- values. 
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Figure 3. Plot of AUeff as a function of ln(L) at 
l3 = 2.4. The curve corresponds to the rescaled 
data of Fig. 2 using the scaling hypothesis Eq. 
10 and Eq. 11. 

5. DISCUSSION 

Based on a scaling hypothesis for a free energy 
difference we have estimated the A/3(/3)-shift in 
the D = 2 0(3)-nonlinear cr-model, see Fig. 4, 
for /3- values as large as /3 = 2.8 on medium sized 
lattices. While at small values of /3 we find agree- 
ment with other numerical determinations, our 




Figure 4. The A/3(/3)-shift (circles) as determi- 
ned by our analysis. The curve corresponds to 
the perturbative result Eq. 3. Circles have been 
connected with a spline. Crosses are from [2]. 

calculation indicates a slower approach to per- 
turbative scaling, than previously acknowledged. 
It would be very helpful, if one could control sca- 
ling relations like Eqs. 10, 11 and possible scaling 
deviations from them, with the help of analytic 
calculations for Ueff(M) in the future. If these 
calculations would confirm our assumptions, then 
possible benefits are significant, as one might then 
more closely approach the perturbative limit of 
asymptotically free theories in numerical simula- 
tions. 
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